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The  "perpetual  salt  fountain"  can  occur  in  an  area  of  the  ocean 
where  the  water  is  vertically  stable.  Warm,  saline  water  overlies  cold, 
relatively  fresh  water.  A  tube  is  placed  in  the  water  from  just  beneath 
the  ocean  surface  to  just  beneath  the  halocline.  Water  is  pumped  up  the 
tube  to  initiate  the  flow.  As  the  water  in  the  tube  is  warmed,  it  be- 
comes less  dense  than  (and  thus  buoyant  with  respect  to)  the  surrounding 
water,  and  the  flow  will  continue  after  pumping  is  stopped. 

The  effect  of  the  vertical  conduction  of  heat  by  turbulent  mixing 
on  the  velocity  of  flow  in  the  "perpetual  salt  fountain"  is  investigated. 
Temperature  and  salinity  profiles  with  depth  from  hydrographic  soundings 
are  considered  to  determine  flow  velocities  in  the  ocean. 

The  calculations  show  that  vertical  heat  conduction  has  a  negligible 
effect  on  the  rate  of  flow  in  the  presence  of  horizontal  heat  conduction, 
the  value  of  the  average  salinity  of  a  column  of  water  minus  the  salinity 
of  the  water  entering  the  bottom  of  the  tube  is  the  primary  factor  in- 
fluencing the  rate  of  flow,  and  the  effect  of  the  temperature -depth  profile 
on  the  velocity  of  flow  is  negligible  so  long  as  initial  stability  re- 
quirements are  met. 
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1.  Introduction. 

Stommel,  Arons,  and  Blanchard  [8]  first  wrote  about  an  "oceano- 
graphic  curiosity"  which  they  called  the  "perpetual  salt  fountain." 
Consider  a  situation  in  which  there  is  a  layer  of  cold,  relatively 
fresh  water  which  underlies  a  warm,  saline  upper  layer  of  water.  The 
density  of  water  is  a  function  of  temperature,  salinity  and  pressure, 
and  the  density  of  the  upper  layer  must  be  less  than  the  density  of 
the  lower  layer  for  stability. 

Now  insert  a  tube  into  the  water  which  extends  from  just  beneath 
the  surface  of  the  water  to  the  top  of  the  lower  layer  of  water.  A 
power  source  is  provided  to  initiate  an  upward  flow  of  water.  As  the 
water  in  the  lower  layer  is  started  up  the  tube,  it  is  warmed  by  heat 
flowing  into  the  tube  from  the  surrounding  warm  water.  Since  the 
water  in  the  tube  is  less  saline  than  the  water  outside  the  tube,  it 
becomes  less  dense  than  the  outside  water  and  is  forced  upward. 

The  flow  process  requires  no  source  of  power  other  than  that  used 
to  initiate  the  flow,  and  flow  will  continue  so  long  as  there  continues 
to  be  a  well-defined  two- layer  system.  The  principle  will  work  for  a 
flow  in  either  direction. 

Groves  [2]  calculated  the  rate  of  this  flow  in  the  tube  to  de- 
termine what,  if  any,  practical  value  this  perpetual  flow  phenomenon 
might  have.  Groves'  calculation  was  based  on  a  simple  model  where  the 
water  temperature  was  considered  a  constant  down  to  the  bottom  of  the 
tube  which  just  extended  into  water  of  a  colder  temperature.  Groves 
determined  that  the  volume  rate  of  flow  for  his  tube  would  be  too  small 
to  be  of  practical  significance  for  fish  farming  (that  is,  use  of  this 

principle  as  a  source  for  providing  nutrient- rich  deep  water  at  the 
surface  to  increase  the  amount  of  edible  fish  in  an  area). 


Groves'  model  neglected  the  effects  of  vertical  heat  conduction 
within  the  tube,  which  should  increase  the  warming  of  the  water  inside 
the  tube  and  provide  a  corresponding  increase  in  rate  of  flow.  If  this 
were  the  case,  then  perhaps  a  larger  tube  could  be  used. 

The  additional  heating  which  is  a  result  of  vertical  heat  con- 
duction would  be  sufficient  to  warm  the  increased  volume  of  water  in 
a  larger  tube.  If  the  flow  velocity  in  a  tube  of  larger  inner  radius  was 
approximately  equal  to  the  flow  velocity  in  the  tube  of  smaller  inner 
radius,  the  tube  with  the  larger  inner  radius  would  provide  an  increase 
in  volume  rate  of  flow. 

It  was  considered  that  a  significant  increase  in  the  velocity  of 
the  water  in  the  tube  might  increase  the  volume  of  flow  sufficiently  to 
warrant  further  investigation  of  this  process  and  its  practical  appli- 
cations. The  effects  of  wall  friction  on  the  temperature  of  the  water 
inside  the  tube  were  considered,  but  the  result  was  so  small  that  the 
effect  was  considered  negligible.  Calculations  of  temperature  change 
due  to  this  stress  are  contained  in  Appendix  C. 

2.  Solution  for  a  model  which  includes  the  effects  of  vertical 
heat  conduction. 

In  order  to  see  the  effect  of  vertical  heat  conduction  on  the 
rate  of  flow  of  the  water  inside  the  tube,  let  us  consider  the  same 
model  that  Groves  worked  with,  but  include  vertical  heat  conduction. 

Assume  that  the  water  density  is  a  linear  function  of  temperature 
and  salinity,  and  that  the  water  is  incompressible.  Therefore,  the  rate 
of  change  of  density  with  temperature,  ?/£0  ,  and  the  rate  of  change 
of  density  with  salinity,  %S  ,  are  constant,  a  valid  assumption 
within  the  small  range  of  temperature  and  salinity  considered.  Then, 


where  ^  is  the  water  density  and  0  and  S  respectively  are  the  temp- 
erature and  salinity  of  the  water.  The  heat  flow  into  the  tube  is 
calculated  by  the  equation 

g  ktWf)  *  <*  2fl*d|  <-£*>*■  «^  d^  (2) 

where  D/Dt  designates  the  material  derivative  following  the  mean 
motion  of  the  water,  q  is  the  heat  per  unit  mass,  n  is  the  heat  con- 
ductivity of  tube  material  (say,  copper  ),  r  is  the  inner  radius  of 
the  tube,  0(y)   is  the  temperature  of  the  water  in  the  tube,  w  is  the 
wall  thickness  of  the  tube,  and  K  is  the  coefficient  of  vertical  heat 
conduction  for  turbulent  flow  calculated  in  Appendix  D.  At  the  sea 
surface  y=0,  and  y  is  positive  downward.  If  we  assume  a  steady  state, 

Dt    *    «ty    <H  "  <ty  (3) 

since  the  local  change      %tsO.   (  The  average  flow  velocity  is  v.) 
Combining  the  constant  terms,  equation  (2)  becomes 

d!$H-$^-o<<p=-o<_e  (4) 

where  £=hcM?/LK  and  cX=2n/Krw.  The  specific  heat  of  sea  water  is  c,  h  is 
the  depth  of  the  bottom  of  the  tube,  and  L  is  the  length  of  the  tube. 

The  general  solution  of  the  homogeneous  equation  associated 
with  (4)  is 

A.y-    .     -r   -  A> 


&(*)-  T,eA*+  Tae>* 


(5) 


where  J^  and  J2  are  undetermined  constant  coefficients, 


and  , 


Az~   J 


h 


(7) 


(9) 


2.1.  Solutions  for  fourth-power  temperature  and  salinity  profiles. 

In  order  to  find  a  particular  solution  to  equation  (4) ,  a  fourth- 
power  temperature  profile  was  assumed: 

©6/)=A-C^.  (8) 

This  function  was  chosen  because  it  can  be  easily  fit  to  actual 
temperature-depth  curves.  The  value  of  sea-surface  temperature  can 
be  inserted  for  A,  and  the  value  (Ts-  T.  )/h  inserted  for  C.  The 
complete  solution  for  equation  (4)  is  __ 

The  water  entering  the  bottom  of  the  tube  is  at  the  same  temper- 
ature as  the  surrounding  water  at  that  depth;  therefore,  0(h)*9(h). 
For  steady  flow,  the  water  at  the  top  of  the  tube  is  assumed  to  be  at 
the  same  temperature  as  the  water  surrounding  it.  This  is  a  reasonable 
assumption  since  we  are  considering  vertical  heat  conduction,  and  the 
water  at  the  top  of  the  tube  undergoes  vertical  turbulent  mixing  with 
the  surrounding  water.  Thus,  the  other  boundary  condition  is 
0(0) =9(0).  Only  upward  flow  is  considered  in  this  paper. 

Applying  the  boundary  conditions  to  equation  (9)  and  solving  for 
J-i  and  Jy  t  — I 

b(t^-))-EU-Fh2-GU3l  eA,f  m 

+■  a-  d-e>- jv-sy-c/  , 

The  constants  D,  E,  F,  and  G  are  functions  of  oCt  $   ,  and  C  for  a  given 
velocity: 
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0 


p  =  is 1 2  c<x.  + 12  e2c/  5 

Now  assume  that  the  water  is  removed  from  the  top  of  the  tube  as 
fast  as  it  flows  up  so  that  no  additional  pressure  head  is  formed* 
Then  the  pressure  difference  driving  the  flow  is  the  hydrostatic  pres- 
sure outside  the  tube  at  y=»h,  less  the  hydrostatic  pressure  in  the  tube 
at  this  depth. 


Equation  (1)  gives  an  internal  pressure  difference  of 
K      nh 


e< 


+  (%\sfo<%l^l 


4        Jo 

and  an  external  pressure  difference  of 


(11) 


*    . 


(12) 


Since  both  pressures  are  atmospheric  at  the  surface,  subtraction  of 
equation  (12)  from  equation  (11)  gives  the  pressure  which  drives  the 


flow:  u 


Ap-| 


s^ysi^fs  ^)-^EJ  ^  0 


(13) 


Dfcnsity  is  constant  to  within  one  per  cent,  and  from  the  continuity 
equation.  ^  v=constant.  This  implies  that  v-constant  and.  in  fact, 
will  be  assumed  constant.     Then  the  wall  friction #Tt  is  also  constant: 


T-if60?v*. 


8 

Here.  f(R)  is  called  the  "resistance  coefficient.11  (Rouse  [5]).     Due 
to  the  small  value  of  viscosity  for  sea  water,  the  flow  will  be  as- 
sumed to  be  turbulent,  and  the  relationship 


(14) 


11 


f(fO  -  .3l<°t/f<l/i 


(15) 


for  smooth  tubes  will  be  used.  Here,  R  is  the  Reynolds  Number, 

R  =  2\r/r/*  (16) 

where  ^  is  the  kinematic  viscosity  of  sea  water. 

Since  the  water  in  the  tube  is  not  accelerating,  the  driving 
force  due  to  pressure  must  equal  the  retarding  force  due  to  wall 
friction: 

2  L  i  (17) 


AT  - 


Combining  equations  (14)  and  (17)  gives 


H? A-p-o. 

Solving  explicitly  for Zrp  in  equation  (13)  necessitates  prescribing 
a  salinity- depth  curve*     A  fourth-power  salinity  profile 

SCi)='Q-Pt*  (19) 

gives  the  salinity  curve  with  depth  the  same  physically  reasonable 
form  as  the  temperature  curve.  The  surface  salinity  is  Q  and 
P=(Sc=SjJ/h  ,  Equation  (13)  becomes  rrt 

if  |(ieh-eFh^gM)/>,0-e^K)-^(i-e^K)f) 

Equation  (20)  is  difficult  to  evaluate  in  its  present  form,  but 
there  are  several  approximations  which  can  be  made  to  simplify' the  ex- 
pression without  sacrificing  accuracy. 
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First,  assume  that  the  rate  of  flow  will  be  less  than  ten  meters 
per  second.  Calculated  values  were  much  less  than  this  value,  and  a 
velocity  of  ten  meters  per  second  would  require  an  extremely  large 
pressure  force  to  drive  the  flow. 

Second,  assume  that  the  tube  thickness  is  very  much  less  than  the 
inner  radius  of  the  tube.  A  thin  tube  is  required  to  allow  the  water 
in  the  tube  to  reach  approximately  the  same  temperature  as  the  water 
outside  the  tube  as  it  flows  upward. 

Third,  assume  terms  that  are  10"  (order  of  magnitude)  as  large 
as  the  dominant  terms  are  negligible. 

From  equation  (5),  Aj ,  is  a  function  of  c<  and  y^  and  will  always 
be  a  negative  number.  If  h  is  greater  than  (say)  50  meters,  (it  will 
usually  be  much  greater  than  this),   /\|h  is  a  large  negative  number, 
and  e  '  is  negligible. 

The  constants  D,  E,  F,  and  G  are  such  small  numbers  that,  even 

2      3 
when  multiplied  by  1,  h,  h  ,  and  h  respectively,  they  are  still  negli- 
gible compared  to  other  terms. 
Equation  (20)  becomes 


and  combining  equations  (18)  and  (21)  gives 


(21) 


(22) 


This  is  a  closed  form  of  the  solution  for  the  flow  velocity  which  in- 
cludes the  effects  of  vertical  heat  conduction. 

If  vertical  heat  conduction  is  not  considered,  the  governing 
differential  equation  (2)  simplifies  to 


13 


£(^d,)--^v^ 


which  has  the  solution 


(23) 


<p60=tfe** 


Pe^Te(r)dT^^h)e'^h-^  (24) 


where  &  =2nL/wcvrho  , 

Solving  equation  (24)  for  the  fourth-power  temperature  profile  gives 

" £      l~r   i*'  — ^~J ' 

Using  the  fourth-power  salinity  profile  and  equation  (25)   for  0(y), 
equation  (13)  becomes 

where  negligibly  small  terms  have  been  ignored.     The  solution  for  fourth- 
power  temperature  and  salinity  profiles  when  vertical  heat  conduction  is 
not  considered  is 

Equations   (22)   and  (27)  were  solved  using  the  temperature  and  sali- 
nity profiles  pictured  in  figure  1,  using  identical  tube  dimensions,, 
The  flow  velocity  which  was  influenced  by  vertical  heat  conduction  was  less 
than  one-tenth  of  one  per  cent  greater  than  the  flow  velocity  calculated 
ignoring  vertical  heat  conduction  (Table  1,)„ 
2o2     Solutions  for  linear  temperature  and  salinity  profiles. 

In  order  to  consider  the  effect  of  different  temperature  and  sali- 
nity profiles,  linear  profiles  were  assumed: 


(27) 
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q(%)   =  A\-  A/y,  (28) 

and 


(29) 


(30) 


s60=  q-p*  . 

The  solution  to  equation  (4)  then  becomes 

where  X-Nfy*. 

Solving  equation  (13)  using  the  temperature  profile  of  equation  (30) 
gives 

where  negligibly  small  terms  have  been  ignored.  The  solution  for  linear 
temperature  and  salinity  profiles,  including  vertical  heat  conduction 
effects  is 

Assuming  linear  temperature  and  salinity  profiles,  equation  (24)  becomes 
and  the  solution  to  equation  (13)  is 


^  X    I  (34) 


where  negligibly  small  terms  have  been  ignored.  Thus,  the  solution  for 
flow  velocity  for  linear  temperature  and  salinity  profiles  where  verti- 
cal heat  conduction  is  ignored  is 
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Equations  (32)  and  (35)  were  solved  using  the  temperature  and  sali- 
nity profiles  pictured  in  Figure  2,  with  identical  tube  dimensions,,  Again, 
the  rate  of  flow  influenced  by  vertical  heat  conduction  was  less  than  one- 
tenth  of  one  percent  greater  than  the  rate  of  flow  calculated  ignoring 
vertical  heat  conduction*  (Table  1.). 

It  is  interesting  to  compare  the  actual  calculations  of  flow  veloci- 
ty for  the  profiles  in  Figure  1  with  those  related  to  Figure  2.  Vertical 
heat  conduction  was  included,  but  this  is  not  pertinent  to  the  discussion 
below.  The  velocity  solutions,  equations  (22),  (25),  (32)  and  (34)  all 
contain  the  term  lQ-Ph/2-S(hJ  which  is  the  average  salinity  of  the  water 
between  the  surface  and  the  depth  y=h,  "s",  less  the  salinity  at  y=h,  S(h) . 
In  making  these  calculations  it  was  noted  that  the  main  factor  causing 
the  one  centimeter  per  second  difference  in  rate  of  flow  between  the  pro- 
files in  Figures  1  and  2  was  the  value  of  the  S-S(h)  term. 

Calculating  the  value  of  S-S(h)  for  the  salinity  profiles  in  Figure 
1  and  in  Figure  2  showed  that  the  two  values  were  almost  identical. 
However,  the  value  of "S* obtained  by  applying  the  fourth-power  salinity 
profile  approximation  to  the  salinity  curve  of  Figure  1  gave  an  erroneous 
value  forl^S(h).  The  error  in  the  calculation  of  S  by  substituting  an 
approximate  function  for  the  actual  salinity  curve  caused  a  marked  error 
in  flow  velocity  result.  This  error  was  later  verified  by  the  computer 
program  determination  of  flow  velocity  which  required  no  approximation  of 
IT.  This  indicates  the  great  importance  of  the  salinity  profile  which 
will  be  discussed  in  more  detail  below. 
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The  flow  rates  for  the  0-S  profiles  contained  in  Figures  1  and  2  were 
also  calculated  utilizing  a  digital  computer  (see  below),  Here  the  ef- 
fects of  vertical  heat  conduction  were  not  included.  The  rate  of  flow  in 
each  case  was  within  one  per  cent  of  the  other,  and  was  essentially  the 
same  as  the  flow  rate  calculated  by  hand  from  equation  (32) . 

TABLE  1 
RESULTS  OF  CALCULATIONS  FOR  RATE  OF  FLOW 


Case 
No. 

6(y)   $  S(y) 
Profiles 

"Fixed  # 
Parameters" 

Hand  ## 
Calculations 

Computer  ## 
Calculations 

1 

Figure 

1 

Approximate 

17.32 

2 

Figure 

2 

Approximate 

17, 

36 

17.30 

3 

Figure 

1 

Approximate 

16, 

,32* 

4 

Figure 

2 

Approximate 

28, 

52 

27.88 

5 

Figure 

1 

Precise 

17.40 

6 

Figure 

3 

Precise 

22.68 

7 

Figure 

4 

Precise 

26.60 

8 

Figure 

5 

Precise 

20.89 

9 

Figure 

6 

Precise 

8.03 

10 

Figure 

5 

Precise*** 

30.94 

11 

Figure 

3 

Precise*** 

33.59 

*  The  flow  velocity  result  is  for  the  hand  calculation  of  the  fourth- 
power  approximation  of  temperature  and  salinity  profiles  illustrated  in 
Figure  1. 

**  The  inner  radius  of  the  tube  was  increased  to  20  centimeters  and  the 
thickness  of  the  tube  increased  to  two  centimeters. 

***  The  tube  was  inserted  vertically  so  that  the  length  of  the  tube  equals 
the  depth  of  the  bottom  of  the  tube, 

#  Appendix  E  contains  "fixed  parameter"  values.  Except  where  otherwise 
noted,  the  inner  radius  of  the  tube  is  0,1  meters  and  the  wall  thickness 
is  0.002  meters. 

##  Flow  velocity  results  are  in  centimeters  per  second. 
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Marked  differences  in  rate  of  flow  which  occur  when  the  tube  is 
inserted  into  water  which  has  different  temperature  and  salinity  pro- 
files are  a  direct  result  of  the  difference  in  the  salinity  profile  from 
one  location  to  another.  The  temperature  profiles  did  not  significantly 
influence  the  rate  of  flow  within  the  small  but  realistic  range  of 
temperature  and  salinity  used  by  the  author. 
2,3  Effects  of  vertical  heat  conduction. 

The  importance  of  vertical  heat  conduction  was  investigated  in 
another  manner.  The  wall  thickness  of  the  tube  was  increased  by  a  factor 
of  ten.  The  flow  velocity  decreased,  but  the  decrease  was  less  than  one 
per  cent.  Increasing  the  inner  radius  of  the  tube  from  ten  centimeters 
to  20  centimeters  and  retaining  the  thicker-walled  tube  resulted  in  a 
significant  increase  in  rate  of  flow  of  37%, 

Comparison  of  the  results  of  the  above  calculations  indicated  that 
the  vertical  heat  conduction  due  to  vertical  turbulent  mixing  was  suf- 
ficient to  warm  the  water  in  the  tube.  In  order  to  determine  the  extent 
of  the  warming  due  to  vertical  heat  conduction,  the  velocity  of  flow  was 
then  determined  for  the  situation  in  which  no  horizontal  heat  conduction 
was  allowed.  In  this  case,  the  temperature  in  the  tube  is  a  function 
only  of  vertical  heat  conduction.  The  governing  differential  equation  (2) 
simplifies  to 

which  has  the  solution 
where  a=KL/cvh^>  , 


18 


Solving  equation  (13)  using  equation  (26)  and  linear  temperature  and 
salinity  profiles, 

where  negligibly  small  terms  have  been  omitted. 

The  resulting  flow  velocity  is  nearly  identical  to  the  flow  velocity 
calculated  using  the  same  tube,  but  allowing  horizontal  heat  conduction. 
The  process  of  heat  conduction  in  the  vertical  by  turbulent  mixing  would 
require  a  longer  period  of  initial  pumping  to  reach  a  steady  state  con- 
dition, but  the  results  indicate  that  the  vertical  heat  conduction  pro- 
cess is  sufficient  to  heat  the  water  inside  the  tube. 

The  initial  tube  placement  was  at  an  angle  of  depression  of  30  degrees 
to  the  water  surface.  This  was  done  so  that  the  water  in  the  tube  would 
require  a  longer  time  interval  to  reach  the  surface  and  undergo  increased 
warming.  Since  the  water  in  the  tube  is  sufficiently  warmed  by  vertical 
heat  conduction  to  be  less  dense  than  the  water  surrounding  it  outside  the 
tube,  the  tube  placement  was  changed  to  a  vertical  position.  The  decrease 
in  tube  length  resulted  in  a  corresponding  decrease  in  the  effect  of  wall 
friction  and  a  481  increase  in  rate  of  flow.  Such  an  increase  is  certainly 
significant  and  will  increase  the  volume  rate  of  flow. 
3.  Computer  solutions. 

In  order  to  determine  the  rate  of  flow  for  a  "perpetual  salt  fountain" 
situated  in  the  ocean,  actual  oceanographic  temperature  and  salinity  data 
for  a  given  area  must  be  known.  A  program  to  obtain  flow  velocity  so- 
lutions which  neglected  vertical  heat  conduction  was  written  for  a  CDC 
1604  digital  computer  using  FORTRAN  60  programming.  Using  this  program, 
flow  velocity  was  easily  determined  for  any  area  in  which  temperature  and 
salinity  data  were  available. 
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Temperature  and  salinity  values  for  each  one-meter  interval  over  the 
entire  length  of  the  tube,  tube  dimensions,  and  "fixed  parameters"  are  en- 
tered on  data  cards  in  the  program.  "Fixed  parameters"  include  kinematic 
viscosity  of  sea  water  (>$ ) ,  density  of  sea  water  (P),  specific  heat  of 
sea  water  (c) ,  acceleration  of  gravity  (g) ,  rate  of  change  in  density  with 
change  in  salinity  (??/&S)   and  rate  of  change  in  density  with  change  in 
temperature  (r?/d&) »  Appendix  E  contains  the  values  used  for  "fixed 

parameters ." 

Subroutine  INTE  computes  the  temperature  profile  in  the  tube,  and 

subroutine  DELTP  computes  the  pressure  force  which  drives  the  flow.  Sub- 
routine INTE  provided  some  interesting  problems  due  to  the  size  of  the 
exponents  of  the  exponential  functions.  A  discussion  of  these  problems, 
and  the  solution  to  them  is  included  with  the  computer  program  in 
Appendix  A.  Subroutine  DELTP  is  a  trapezoidal  integration  which  is  as 
accurate  as  can  be  obtained  from  the  small  number  of  data  points  gene- 
rally given  in  a  hydrographic  sounding , 

Initial  data  used  in  the  program  was  set  up  to  treat  the  same  case 

as  was  calculated  by  Groves*  Temperature  and  salinity  profiles  were 

devised  so  that  values  of  9(0) -0(h),  *S-S(h),  and  the  "fixed  parameters" 

would  be  identical  to  those  used  by  Groves.  Additional  test  cases  were 

run  to  compare  with  the  results  of  hand  calculations. 

Values  for  the  constant  terms ty& and  %s  were  then  calculated  from 

the  equation  of  state  given  by  Eckart  [1]  (see  Appendix  B)  and  more  pre- 
cise values  for  the  other  "fixed  parameters"  were  obtained  and  used  with 
data  from  actual  hydrographic  soundings.  Several  cases  of  actual  data 
were  employed  to  get  an  idea  of  the  range  of  values  of  flow  velocity  which 
could  be  expected  in  the  ocean.  The  results  are  contained  in  Table  1,  and 
the  temperature -depth  and  salinity-depth  curves  are  pictured  in  Figures  1-6 , 
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In  the  actual  ocean  data  cases,  the  depth  to  which  the  tube  was  ex- 
tended was  determined  by  the  depth  of  the  halocline.  The  tube  was  always 
extended  into  the  water  just  below  a  negative  halocline, 
4.  Results, 
4,1  Vertical  heat  conduction  versus  no  vertical  heat  conduction. 

Initial  flow  velocity  calculations  were  made  by  hand  using  the  initial 
values  of  the  "fixed  parameters"  contained  in  Appendix  E.  The  rate  of 
flow  for  identical  tube  dimensions  was  calculated  for  0-S  profiles  con- 
tained in  Figures  1  and  2, 

The  calculation  of  flow  rate  was  made  twice  for  the  linear  9-S  pro- 
files of  Figure  2,  once  including  the  effects  of  vertical  heat  conduction 
and  once  neglecting  the  effects  of  vertical  heat  conduction.  In  each 
case,  inclusion  of  vertical  heat  conduction  had  a  negligible  effect  on 
the  rate  of  flow. 

Flow  velocity  calculations  were  then  made  for  the  0-S  profiles  of 
Figure  1  using  the  fourth-power  temperature  and  salinity  profiles.  These 
functions  were  selected  because  they  approximated  the  physical  appearance 
of  the  curves  in  Figure  1  and  because  they  were  easily  integrable.  Flow 
rate  solutions  for  the  vertical  heat  conduction  and  no  vertical  heat  con- 
duction equations  were  again  identical  to  four  significant  figures. 

Flow  velocity  was  then  calculated  for  a  situation  in  which  no  hori- 
zontal conduction  of  heat  was  allowed,  so  that  only  vertical  heat  con- 
duction was  possible.  The  rate  of  flow  for  this  calculation  was  nearly 
identical  to  the  rate  of  flow  calculated  for  two  other  cases  which  used 
the  same  tube  dimensions,  "fixed  parameter"  values,  and  0-S  profiles. 
The  first  allowed  both  horizontal  and  vertical  heat  conduction,  and  the 
second  allowed  only  horizontal  heat  conduction. 
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On  the  basis  of  these  calculations,  it  was  concluded  that  the  effect 
of  vertical  heat  conduction  on  the  rate  of  flow  of  the  water  inside  the 
tube  was  negligible  in  the  presence  of  horizontal  conduction ,  Thus  it 
was  decided  to  make  use  of  the  computer  program  which  did  not  include 
the  effects  of  vertical  heat  conduction  for  the  remaining  calculations. 
4.2  Computer  results  versus  Groves'  result. 

Utilizing  the  computer  program  allowed  exact  values  of  temperature 
and  salinity  at  each  depth  to  be  used.  Temperature  and  salinity  values 
from  Figures  1  and  2  were  used  in  the  program  with  a  resulting  rate  of 
flow  of  17* 30  centimeters  per  second  (cps)  for  Figure  1  profiles  and 
17.32  cps  for  Figure  2  profiles.  The  17.32  cps  flow  velocity  associated 
with  the  profiles  of  Figure  2  was  in  close  agreement  with  the  rate  of 
flow  calculated  by  hand  for  the  same  profiles.  The  17.30  cps  rate  of 
flow  calculated  for  the  profiles  in  Figure  1  was  much  greater  than  the 
rate  of  flow  calculated  by  hand  for  the  same  profiles.  Investigation  of 
the  hand  calculation  revealed  that  the  S-S(h)  term  calculated  by  hand 
differed  significantly  from  the  actual  value  of  S-S(h)  obtained  from 
Figure  1.  This  difference  was  attributed  to  use  of  the  fourth-power 
salinity  profile  as  an  approximation  of  the  actual  salinity  profile 
when  making  hand  calculations.  S-S(h)  calculated  from  the  curve  of  the 
fourth-power  salinity  profile  is  significantly  different  from  the  true 
value  of  S-S(h)  in  Figure  1.  The  influence  of  the  approximation  to  the 
form  of  the  temperature  profile,  however,  was  determined  to  be  negligible. 

The  profiles  of  temperature  and  salinity  in  Figures  1  and  2  were 
constructed  so  that  the  values  of  S-S(h)  and  9(0) -9(h)  would  be  identical 
to  those  selected  by  Groves.  Tube  dimensions  and  "fixed  parameter" 
values  were  also  those  used  by  Groves.  Groves'  calculated  rate  of  flow 
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was  17.50  cps  as  compared  with  a  rate  of  flow  of  17.32  cps  obtained  from 
the  computer  program.  The  difference  in  flow  velocities  can  be  attribu- 
ted to  a  rounding-off  procedure  used  by  Groves  to  simplify  calculations. 

The  values  of  %a  and  ^S  calculated  in  Appendix  B  were  then 
inserted  into  the  original  program  along  with  more  precise  values  of  the 
other  "fixed  parameters"  to  determine  what  effects  the  more  precise  values 
would  have  on  the  rate  of  flow.  The  resulting  velocity  increase  of  eight 
centimeters  per  second  was  not  significantly  greater,  but  the  more  pre- 
cise values  for  the  'fixed  parameters"  were  retained  for  the  remaining 
temperature  and  salinity  profiles. 

Temperature  and  salinity  data  from  actual  hydrographic  soundings  were 
also  used.  Figures  3  through  6  illustrate  the  actual  profiles  used.  It 
was  verified  that  areas  of  most  pronounced  negative  haloclines  produced  the 
greatest  rates  of  flow. 

A  marked  increase  in  rate  of  flow  was  noted  when  the  tube  was  in- 
serted vertically  instead  of  at  a  depression  angle  of  30  degrees  to  the 
ocean  surface.  The  481  increase  in  speed  of  flow  will  provide  an  equal 
increase  in  volume  rate  of  flow. 

Rate  of  flow  was  also  increased  by  increasing  the  inner  radius  of  the 
tube.  This  is  doubly  significant  since  an  increase  in  flow  velocity  in- 
creases volume  rate  of  flow,  and  increased  tube  radius  will  also  increase 
tlie  volume  rate  of  flow, 
5.  Conclusions  and  acknowledgements. 

The  test  cases  discussed  in  this  paper  have  demonstrated  some  interest- 
ing features.  First,  the  temperature -depth  structure  is  of  relatively 
little  importance.  The  surface  layer  must  be  sufficiently  warm  to  in- 
sure stability,  but  the  magnitude  of  the  temperature  change  between  the 
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surface  and  the  bottom  of  the  tube  has  a  negligible  effect  on  the  velo- 
city of  flow. 

Second,  the  primary  parameter  influencing  the  rate  of  flow  of  water 
in  the  tube  is  the  difference  between  the  average  salinity  in  the  layer 
through  which  the  tube  extends,  and  the  salinity  of  the  water  entering 
the  bottom  of  the  tube0  To  calculate  the  rate  of  flow  for  a  given  case 
to  within  one  per  cent  of  its  precise  value,  a  salinity- depth  curve  must 
be  known,  A  simple  are a -averaging  process  can  be  used  to  determine  Sp 
the  average  salinity,  and  S(h)  is  known.  Then, 

Third,  equation  (39)  gives  excellent  results  because  the  effects 
of  vertical  heat  conduction  heat  the  water  in  the  tube  sufficiently  to 
make  it  buoyant ,  Although  vertical  heat  conduction  does  not  directly 
increase  the  velocity  of  flow,  the  warming  effects  allow  the  substi- 
tution of  a  poorer  heat-conducting  material  for  the  copper  tube  A 
metal  which  is  resistant  to  sea  water  corrosion  would  be  a  better  material 
for  the  tube  In  addition,  a  relatively  thick  tube  with  a  larger  inner 
radius  may  be  used.  This  will  increase  the  volume  rate  of  flow  and 
provide  a  more  durable  tube. 

For  example,  using  data  from  the  equatorial  Atlantic  Ocean  region 
in  case  #11  (Figure  3),  the  volume  rate  of  flow  is  10 , 5  liters  per  second,, 
This  is  a  901  increase  in  volume  rate  of  flow  over  that  calculated  by 
Groves,  When  the  inner  radius  of  the  tube  was  enlarged  in  case  #4 
(Figure  2),  the  volume  rate  of  flow  increased  to  34,6  liters  per  second,, 
This  is  an  increase  of  5301,  which  is  definitely  significant. 

It  is  therefore  concluded  that  the  "perpetual  salt  fountain"  is  more 
than  an  "oceanographic  curiosity,"  Upward  flow  of  nutrient-rich  water 
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from  well  below  the  surface  in  sufficient  volume  could  make  a  fertile 
fishing  area  out  of  an  area  which  is,  at  present,  barren.  A  relatively 
inexpensive  and  simple  method  for  providing  this  nutrient-rich  water, 
such  as  the  method  described  in  this  paper,  is  worthy  of  additional 
experimentation . 

The  author  gratefully  acknowledges  the  guidance,  encouraging 
assistance  and  critical  comments  of  Assistant  Professor  T.  Green.  My 
thanks  also  go  to  Mr.  R.  D.  Brunell  of  the  United  States  Naval  Post- 
graduate School's  computer  facility  staff  for  his  valuable  programming 
assistance. 
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APPENDIX  A 

/ 
COMPUTER  PROGRAM  FOR  FLOW  VELOCITY  CALCULATIONS 


JOR0532FM,HINMAN,K.G. ♦FLOWUP 
PROGRAM    FLOWUP 

DIMENSION   T(401) ,    5(401),    PHI(401) 

COMMON  NVT,H,GAMMA,DELH,PHI,S,T,SALTH,THETH,G,DELTP»CHDS,CHDT 
READ  50,  NVT 
50  FORMAT  (1215) 

READ  10»  (T(  I  )  ,S( I ) , I  =  1»NVT) 

10  FORMAT  (6E12.6) 

READ  10»  R»E, CON, PL, W,D,SH, DEPTH, V0»EP1,THETH 
READ  10,  H,G,DELH,SALTH,CHDS,CHDT 
PRINT  12,  R,E,CON,PL,W,D 
PRINT  12,  SH,DEPTH,V0,EP1»THETH 
PRINT  12,  H,G,DELH,SALTH,CHDS,CHDT 
12  FORMAT  (6E12.6) 
V  =  V 

19  V=V+0.01 
RE=2.0*V*R/E 

GAMMA=4.0*CON*PL/(W*D*SH*E*RE*DEPTH) 
CALL  INTE 

CALL  FINAL 
IF(DELTP-O.O)  55,55,56 

55  V=V-0.01 
GO  TO  57 

56  TAU=DELTP*R/(2.0*PL) 
IF(RE-2000.0)  20,20,21 

20  FUNC=64.0/RE 
GO  TO  22 

21  FUNC=0.3164/( (RE)**( 1.0/4.0) ) 

22  VFL=(8.0*TAU/(FUNC*D) )**( 1.0/2.0)  I 

11  FORMAT  (30H   THE  UPWARD  FLOW  VELOCITY  IS  2E20.6) 
IF(VEL-V)  55,24,19 

57  V=V+0.001 
RE=2.0*V*R/E 

GAMMA=4.0*CON*PL/(W*D*SH*E*RE*DEPTH) 
CALL  INTE 

CALL  FINAL 
IF(DELTP-O.O)  40,40,41 

40  V=V-0.001 
GO  TO  42 

41  TAU=DELTP*R/(2.0*PL) 
IF(RE-2000.0)  43,43,44 

43  FUNC=64.0/RE 
GO  TO  45 

44  FUNC  =  0.3164/(  (RE)**(0.25)  ) 

45  VEL=(8.0*TAU/(FUNC*D> )**(0.5) 
IF(VEL-V)  40,24,57 

42  V»V+0.0001 
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RE=2.0*V*R/E  / 

GAMMA=4.0*CON*PL/(W*D*SH*E*RE*DEPTH) 

CALL  INTE 

CALL  FINAL 

IF(DELTP-O.O)  46,46,47 

46  V=V-0.0001 
GO  TO  7a 

47  TAU=DEL7P*R/(2.0*PL) 
IF(RE-2Q00.O)  71,71,72 

71  FUNC=64.0/RE 
GO  TO  73 

72  FUNC  =  0.3164/(  (RE)**(  1.0/4.0).)  f 

73  VEL=(8.0*TAU/(FUNC*D) )**( 1.0/2.0)    . 
IF(VEL-V)  46,24,42 

70  V=V+0. 00001 
RE=2.0*V*R/E 
PRINT  60,  RE 

60  FORMAT  (25H   THE  REYNOLDS-  NUMBER  IS  E20.6) 
GAMMA=4.0*CON*PL/(W*D*SH*E*RE*DEPTH) 
PRINT  61,  GAMMA 

61  FORMAT  (14H  GAMMA  EQUALS  E20.6) 
CALL  INTE 

CALL  FINAL 
PRINT  63.  DELTP 
63  FORMAT  ( 14H  DELTP  EQUALS  E20.6) 
TAU  =  DELTP*R/ ( 2.0*PL  ) 
PRINT  62»  TAU 

62  FORMAT  (12H  TAU  EQUALS  E20.6) 
IF(RE-2000.0)  74,74,75 

74  FUNC=64.0/RE 
GO  TO  76 

75  FUNC  =  0.3164/( (RE)**( 1.0/4.0)  ) 

76  VEL=(8.0*TAU/(FUNC*D) )**( 1.0/2.0) 
PRINT  11,  VEL,V 

DIFF=ABSF( (VEL-V)/V) 

IF(DIFF-EPl)  24,24,70 
24  STOP 

END 

SUBROUTINE  INTE 

DIMENSION  T(401),  S(401)»  PHK401) 

COMMON  NVT, H, GAMMA, DELH,PH I, S,T»SALTH»THETH,G, DELTP, CHDS»CHDT 

DO  10  J=1,NVT 

PHI (J)=0.0 
10  CONTINUE 

IJ=NVT 

PHI ( IJ)=THETH 

BIG=GAMMA+1.0 

BIGEX=EXPF(GAMMA) 

BNE=GAMMA-1.0 
20  IJ=IJ-1 

K=IJ+1 

XI=(T(IJ)-(T(K)*BIG) )/(GAMMA*BIGEX) 
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XJ=( (T( IJ)*BNE)+T(K) ) /GAMMA 

XK=PHI (K) /BIGEX 

PHI  (  IJ)=,XI+XJ+XK. 

IF(IJ-l)  30,30*20 
30  RETURN 

END 

SUBROUTINE  FINAL 

DIMENSION  T(401),  S(401),  PHI (401) 

COMMON  NVT.H»GAMMA.DELH,PHI»S,TtSALTH.THETH»G»DELTP,CHDS,CHDT 

TOTAL=0.0 

1  =  1 

PRESS=(S( I)-SALTH)*CHDS+(T(I J-PHH I ) ) *CHDT 
40    1=1+1 

SUMP=( (S( I )-SALTH)*CHDS)+< (T( I )-PHI ( I ) ) *CHDT ) 

TOM=DELH*(PRESS+0. 5* (SUMP-PRESS) ) 
PRESS=SUMP 
T0TAL=T0TAL+T0M 
IF(NVT-I)  50,50*40 
50  DELTP=G*TOTAL 
RETURN 
END 
END 


~> 
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The  computer  program  is  an  iteration  process  for  obtaining  the 
velocity  of  the  water  in  the  tube.  It  was  written  to  solve  actual 
oceanographic  cases  using  the  form  of  the  solution  given  by  Groves  [2] , 

Subroutine  DELTP  is  utilized  to  compute  the  pressure  integral, 

h 

The  subroutine  is  a  trapezoidal  integration  which  assumes  a  linear 
relationship  between  successive  points.  No  greater  accuracy  is  desirable 
because  hydrographic  soundings  contain  salinity  and  temperature  data  for 
relatively  few  depths,  and  the  linear  relationship  is  generally  assumed 
when  plotting  the  original  data. 

Subroutine  INTE  is  designed  to  compute  the  temperature  of  the  water 
at  each  one-meter  interval  of  depth  inside  the  tube  using 

4>6|)  =  3e**  \e(r)  e-*TdT  +  eWft""^.  <"> 
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Equation  (A2)  relates  water  temperature  in  the  tube  to  depth,  but  it 
cannot  be  solved  in  this  form  on  the  computer  since  the  capacity  of  the 
CDC  1604  is  exceeded  by  exponential  functions  with  exponents  of  magni- 
tude greater  than  709.  A  different  form  of  equation  (A2)  was  necessary 
for  the  computer  program.   -A 
Let 


%(%)*  tie*** Le^T)e~    C'T  (A3) 

(A4) 


^     cp60=  26)+  YM-  CA5) 

In  order  to  get  l(y-Ah) ,  we  let  u, 

u 


In  order  to  get  Z(y-/\h),  we  let  u, 


(A6) 
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and  then  from  equations  (A3)  and  (A6) 

*6-aM=  x6r*^  +  £~*AK  *<*)•       (A7} 

We  now  assume  that  9(T),  where  T  is  a  "dummy"  variable  of  integration, 
varies  linearly  between  consecutive  one-meter  intervals,  and  limit  the 
interval,  Z\h,  to  one  meter. 

Combining  equations   (A5) ,   (A7) ,  and  (A8) 

Since  0(h)=9(h)   is  a  boundary  condition,  the  first  0(y-Zlh)will  be  cal- 
culated for  a  depth  of  (h-1)  meters.     This  integration  subroutine  will 
function  so  long  as  8<  709. 
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APPENDIX  B 
COMPUTATIONS  OF  ty<&MT>^%S 
The  values  for  the  rate  of  change  of  density  with  changing  tempera- 
ture and  salinity  were  assumed  constant  in  the  solution  of  the  pressure 
integral.  Additional  investigation  of  the  value  of  these  constants  was 
considered  important* 

Assuming  that  the  effects  of  pressure  are  negligible, 

Following  Eckart  [1] , 

(?  +  v.)(Ab-A)  =  K 

where 


(B2) 


A  =  /o'(i 


CB3) 


(B4) 


?    / 

and 

-po  =  S8lO  +  3$T  +0.31-&T*  +■  3  5 

ft  -  rn7IS+OJl2ST-a000Vl&Tx  -(o,03W^O,OOOIt)s1x  /£>6(b5) 

A0  =  constant  =3020.  (B6) 

Combining  equations  (B2)  and  (B3)  gives 


^{/-[^(Ao-^r 


(B7) 


Taking  the  partial  derivatives  of  o  with  respect  to  0  and  S  in 
equation  (B7) ,  and  evaluating  the  result  at  an  average  temperature , 
9=15  c»  an  average  salinity,  S=34°/oo,  and  a  zero  pressure,  gives 
^?/6s  =0.77667kg/m3«o/ooand    <^£=0.19987kg/m3'°C. 
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APPENDIX  C 
WALL  FRICTION  CALCULATION 

It  was  considered  that  the  effects  of  the  stress  on  the  water  by 
the  wall  of  the  tube  might  cause  significant  heating  of  the  water.  The 
work  done  on  the  water  per  unit  mass  per  time  is 


^TT^dy        p/T         KWs^J.  (ci) 


joules  A 
Kg °sec  /. 


The  work  per  unit  mass  is 


ZT  V"         Ll  _       ZTh  /joules  \  (C2) 

Dividing  by  the  specific  heat,  c,  and  converting  joules  to  calories  using 
0.2389  calories /joule  at  15°C, 


ZTL, 


(0.2391)   =  AT      (°C). 


(C3) 


Applying  the  data  used  in  the  solution  of  case  #1,  the  temperature  change 
of  the  water  in  the  tube  due  to  wall  friction  is  0.00025°C  for  a  velocity 
of  17.32  cps  and  a  pipe  length  of  600  meters.  This  change  in  temperature 
is  negligible. 
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APPENDIX  D 

COEFFICIENT  OF  VERTICAL  HEAT 

CONDUCTION  FOR  TURBULENT  FLOW 

The  value  for  the  coefficient  of  vertical  heat  conduction  for  turbulent 
flow,  K,  was  investigated.  From  Schlichting,  [6],  the  coefficient  of  eddy 
conductivity  of  heat,  A  ,  is  at  least  as  large  as  the  coefficient  of  eddy 
viscosity,  Ar,  over  the  entire  pipe*  We  assume  the  turbulence  to  be  ap- 
proximately isotropic,  so  that  turbulent  transfer  along  the  pipe  is  the 
same  order  of  magnitude  as  that  in  any  cross  section. 

A  good  estimate  of  the  effective  value  of  Aq  can  then  be  obtained  by 
calculating  an  average  of  A^-over  a  cross  section.  This  is  easily  done 
using  mixing  length  theory.  The  result  is 


—   _   A"       ± 


n* 


Af  -   Ar  - 


j 


«3 


where  r  is  the  interior  radius  of  the  tube,  z  is  measured  inward  from  the 
wall,  u*  is  the  "shear  velocity"  (see  Schlichting) ,  and  k=0.4. 
Then, 

K  =     Xt    t?  .  (D2) 

This  has  a  typical  value  of  5.0  cal./m*  C*sec.  for  the  rates  of  flow 
and  tube  size  considered. 
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APPENDIX  E 
"FIXED  PARAMETER"  VALUES 

The  expression  "fixed  parameters"  refers  to  values  for  kinematic 
viscosity  of  sea  water  (>}),  conductivity  of  tube  material  (  n) ,  densi- 
ty of  sea  water  (£>)»  specific  heat  of  sea  water  (  c  ),  acceleration  of 
gravity  (  g  ),  rate  of  change  of  density  with  changing  salinity  (^P/dS) 
and  rate  of  change  of  density  with  changing  temperature  (<*?/<}&) .  Two 
sets  of  values  for  the  "fixed  parameters"  have  been  used.  One  set  con- 
tains rounded- off  values  for  hand  calculations,  and  the  other  set  contains 
more  precise  values  for  the  computer  program.  One  could  then  determine  the 
effect  on  rate  of  flow  of  changes  in  the  values  of  the  "fixed  parameters." 
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INITIAL  TEMPERATURE-SALINITY  PROFILES 
CONSTRUCTED  TO  APPROXIMATE  GROVES'  DATA 
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Figure   1. 
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LINEAR  TEMPERATURE-SALINITY  PROFILES 
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